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I. INTRODUCTION 


The free electron laser (FEL) consists of a relativistic electron beam which 
passes through the transverse periodic magnetic field of an undulator. This results 
in the exchange of energy between the electrons and a copropagating optical wave. 
Traveling through the undulator the electrons are acted upon by Lorentz forces 
inducing transverse oscillations. This acceleration causes the electron to emit 
electromagnetic radiation. The wavelength of this radiation depends on the 
electron beam energy thereby providing a tunable light source. An FEL configured 
as an oscillator stores the light in an optical cavity. However, in an FEL 
configured as an amplifier, light 1s created and amplified during a single pass of 
the electrons through the undulator eliminating the need for optics. The FEL has 
provided high efficiency and power, but its one overriding attribute which enhances 
its future in science, industry and defense ts its ability to provide, on demand, 
wavelength tunability. As a tunable source of laser light, the FEL is of interest in 
the infrared and ultraviolet spectrums in the study of solid state physics and 
molecular spectroscopy. Medicine is another region in which the FEL may play 
an important role. In the X ray regime, the FEL may lead to holographic imaging 
of DNA molecules and imaging of live biological microstructures in their normal 
physiological environments. Finally, the future of the FEL in missile defense may 
revolutionize naval shipboard defense systems. 


l 


As DEFICIENCIES OF CURRENT SHIP DEFENSE SYSTEMS 

The current systems which comprise the battle group defense in depth has 
a superior long range capability. However, as the inevitable leaker finds its way 
into the inner sanctum of the battle group and a single ship 1s left to defend itself 
against this highly maneuverable and lethal threat, current capabilities are in need 
of much improvement. If we examine the system which 1s the ships final course 
of defense, the phalanx gun system 1s the weak link in the chain. Analyzing the 


phalanx gun through computer simulations, we can first examine a distribution of 


Based on a dispersion of 0.001 radians 





Figure 1-1. Probability of a hit vs. range for a very accurate phalanx gun 
barrel. 


bullets at ranges from 100 meters to 2000 meters using random angular 


dispersions. Figure 1-1 shows the result for the most accurate barrel system with 
a dispersion of 0.001 radians. The graph shows an asymptotic decay in the 
probability of hitting the target as a function of increasing range. From Figure 1-1 
we find that we can achieve a maximum probability of a hit of only 60% inside 
100 meters against a 25 cm radius target. Further, it is accepted in simulations of 
the phalanx that 5 hits are required to kill an inbound missile. Therefore, the 
probability of a kill becomes 8%. Figure 1-2 is another simulation which 1s based 


ona missile kill whose debris is then tracked to impact. We find that a kill at 800 


** Debris Simulation ** 


V5=1000m/s V,=77 m/s R=810.9m 


30|\side view 


x (m) 


T=-1.7 s 


Vv g=258 m/s 





Figure 1-2. Simulation of a missile fragment from 800 meters tracked to 
a waterline impact. 


meters allows a fragment to reach the ship at the waterline with a kinetic energy 
of 10° joules breaching hull integrity. This range of 800 meters we call the 
"critical range". If we return now to Figure I-l and compare with this critical 
range, we find the best case scenario provides a 5% probability of hitting the threat 
and therefore a 0% probability of a kill prior to it achieving a range that will inflict 
damage on the ship. If we use 20 fragments and randomize their size, trajectory 


and velocity, we find in Figure 1-3 thc: at the critical range of 800 meters 60% of 





Range (m) 


Figure 1-3. Results of tracking 20 missile fragments through various 
ranges to impact. 


the fragments reach the ship. Inside of 500 meters 100% of the fragments hit the 
ship. So, what is the answer for close-in point defense where offensive missiles 
have been maximized to defeat conventional defensive systems? Speed of light 


weapons which can reach farther than the 2 km limit of the phalanx. 


B. LASER CLOSE-IN WEAPONS SYSTEM 

In order to destroy an incoming missile, the laser power required on the 
missile surface is about 2 kW/cm’ over a 5 cm radius spot for a one second 
duration. The power required at the missile 1s then 160 kW. Due to the aerosols 
at sea level, the power from the ship to the target diminishes exponentially by, 
P = P, e~, where P, is the power at the source, & = 5 x 10° km” is the extinction 
coefficient and z 1s the range from the source to the target. Based on a range of 
20 km for this new system, the required power leaving the ship 1s 427 kW and 
equates to a beam director radius of 0.5 m. A tunable FEL utilizing a beam energy 
ymc’ = 50 MeV and peak current I = 300 A can have sufficient average current 
to provide 374 kW at the target. This 1s more than double that required to kill the 
missile. An advantage of a speed of light weapon is the fact that based on the 
missile velocity of v, = 1000 m/s the laser beam takes t = 6.67 x 10” s to reach 
the target. In that amount of time the missile only travels 6.7 cm. Even if the 


missile is maneuvering, it will not out-maneuver a speed-of-light weapon. The 


parameters required to achieve this peak power are achievable, and have been 
demonstrated in existing FELs. But, high average power has not been 
demonstrated. The speed-of-light weapon has the advantage of killing a missile 
early enough to allow a greater numbers of engagements. It also allows the 
Operator to assess his success in defeating the missile in seconds allowing the 
tactician to move on to the next threat. 

Why use an FEL instead of current proven conventional lasers? 
Conventional lasers of the required power utilize gases which are both highly 
corrosive and deadly. Further, with each firing of the weapon you deplete your 
resources, and therefore have a limited magazine. The FEL requires no gases and 
as long as there is power to the ship your magazine is never empty. In a typical 
war-at-sea scenario, the enemy 1s required to overwhelm the defender with large 
numbers of threats to assure a 100% probability of destruction of his target. 
Against a speed-of-light defense coupled with current technology long range anti- 
missile missiles, the number of offensive missiles to assure a kill becomes 


prohibitive. 


Il. FREE ELECTRON LASER THEORY 


A. GENERAL 


Typical lasers essentially take some form of energy to excite a lasing 
medium. During the decay from these excited states, light 1s radiated amplifying 
light of the same wavelength, phase and propagating in the same direction as the 
stimulating energy. From the initial energy, a coherent, monochromatic light is 
produced. By contrast, the FEL produces coherent, monochromatic light from the 
acceleration of relativistic electrons. A magnetic field exerting a force on a 
relativistic electron bends its trajectory and electromagnetic theory then demands 
that the electron emit radiation. For a relativistic electron, the power radiated is 
P = (2/3)(y'e*c/R*) [2], where, e is the charge of the particle, c is the velocity of 
light, R is the bending radius, and y is the energy E/mc*. A relativistic particle 
when accelerated radiates into a narrow cone pointing in the direction of its 
motion. If the electrons are subjected to an alternating magnetic field causing them 
to wiggle or undulate, they will then radiate light along a single axis. As this 
radiation builds up, the electrons radiate in the presence of radiation. Figure 2-1 
illustrates the dynamics by which the electron beam energy is converted into 
coherent radiation. A single electron traversing one period of an undulator 1s 


followed while studying its interaction with the light wave. In frame (a) the short 


wavelength optical electric | d at the electrons’ position 1s a maximum and is 
oriented in the same direction as the long wavelength undulator deflection. This 
orientation decelerates the electron and causes it to give energy to the optical field. 
In frame (b) the electron has traveled one-fourth of an undulator period at a speed 
less than c, allowing the optical wave to advance on the electron by one-fourth of 
an optical wavelength. Here the optical electrical field and the undulator field are 
zero. Except for the force of the undulator, the electron feels no additional force. 
Again in frame (c) the optical electric field decelerates the electron. Finally in 
frame (d) the electron has traversed a full undulator period while the optical wave 
has moved a full optical wavelength over the electron. This idealized phase 
relationship demonstrates the electron and optical interaction at resonance. The 
more complex case where a beam of electrons and the optical wave are in a large, 
number of random phase combinations means that while some electrons are 
decelerated, others are accelerated. A consequence of this complex energy 
exchange is that the electrons form electron bunches which collectively exchange 


energy to the radiation field. 


(a) 


(c) 


Figure 2-1. 





(d) 


Illustration of the electron and optical wave interaction 
beginning with frame (a) (upper left), advance of the optical 
wave through a quarter of the undulator period in frame (b) 
(upper right). Frame (c) (lower left) is halfway through the 
undulator and finally frame (d) (lower right) is full passage 
through one period of the undulator. 


B. BASIC PRINCIPLES. 
We examine a linearly-polarized undulator producing coherent linearly- 


polarized light. The undulator field 1s described by [3] 


B = B (0, sin(k,z), 0] oa 
along the undulator (z) axis where B is the undulator peak magnetic field and 
ky=2t/h, is the undulator wavenumber. The undulator length is L=NA,, where A, 
is the undulator wavelength and MN is the number of undulator periods. The 
undulator provides a periodic transverse deflection to the electrons causing them 
to radiate spontaneously. As this spontaneous emission grows, an electron with 
energy y(‘)mc* couples with the co-propagating radiation and exchanges energy 
with the optical pulse. The optical wave has an electric field with strength £ and 
carrier frequency @ = kc = 2nc/X where k 1s the optical wavenumber, A is the 
optical wavelength, and c is the speed of light. The optical mode created by an 


electron traveling through the undulator can be described by [3] 


| 


= E(z,t) [cos y, 0, 0] ~ 


-7x E = E(z,t) [0, cos yw, 0] 


Op | 


where E(z,t) is the wave amplitude, y = kz - wt + $(z,f), and the optical phase 1s 


represented by $(z,f). These fields can be derived from the vector potential 
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A(x,t) = (E(z,0/k) [sin w, 0, 0]. If the amplitude E and phase 9 are held constant, 


then these fields describe a plane wave traveling in the z-direction. 


C. ELECTRON DYNAMICS AND THE PENDULUM EQUATION 


The Lorentz force acting on an electron traveling through the undulator 1s 





yb) 6 ie ee 2-3 
AF <6 B x (B B.)), 
ON a emai 2-4 
7 =<? Jie 
a= BeBe oe 


where, B. is the undulator magnetic fields described by 2-1, E and B are the 
optical electric and magnetic field described by 2-2. Substitution of these fields 


into 2-3 yields 





: ele - B.) (cos w, 0, 0) + BB(-sin (kz), 0, 0)] . 2-6 
at mc : ° 


Assuming the electrons are relativistic, B, ~ 1, we then have 


aYB) . -£ BB (-sin (kz), 0, 0). 2-7 
dt mc 7 


1] 


Now integrate to obtain the transverse elect’ n velocity 


2a 
B= VK co (kez) Om)”, 2-8 


where we have assumed perfect electron injection, so the constant of integration 
is zero. The undulator parameter K=eB,,.A)/2mmc* is a dimensionless quantity 
characterizing the properties of the undulator, and B_, is the rms field strength in 


each period. Substitution into 2-5 and averaging over an undulator period, gives 





aR lege 2-9 


Taking the time derivative, we get 





, 1+K “)y ; ° 
pp, - CY = 3 -pp 7 _. 2-10 
3 leek 


Substituting the electric field 2-2 into 2-4 and using the value of B in 2-8, we 


have 


or = _— cos(y + kz) = aia cos(G + 9) 2-11 
at yay me y2yme 


where the electron phase is € = (k+k,)z-@t. For phases such that -m/2 < C+o < 


1/2, the electron gains energy from the radiation field, whereas for phases such that 


m/2 < C+o < 3n/2, the electron gives up energy to the radiation field. If we 
differentiate the electron phase twice we get d°C/dt° = c(k+k,) dB/dt. Substituting 


2-10 into 2-11 yields the electron pendulum equation 


d?¢  2ek KE 
ae yom 








cos(G + 6) . 2-12 


which governs the electron motion in the presence of the undulator and radiated 
fields. Introducing the dimensionless time t = ct/L and averaging the motion over 


each undulator period [3], we have 





X= = lal cos ( + 4), shia 
dt 


where the dimensionless’ optical field strength is given’ by 
la| = 4nNeEK(J,(é)-J(E))L/y2 y “mc~, J, and J, are Bessel functions of the first 


kind and €=K’/2(1+K’). We further define the dimensionless phase velocity by 
V(t) = = = L[(k+k,)B. - k) . 2-14 
= : 


When v = 0, the electron experiences resonant undulator and radiation field forces. 


The resonance condition can then be expressed by 4 = A,(1+K*)/2y*. This equation 
shows that the desired operational wavelength may be achieved by altering the 
electron beam energy y, the undulator wavelength A, or the undulator parameter 
K. The values of 4, and K are limited by undulator technology; hence most FELs 
currently operate in the infrared or visible regime. However, with the large value 
of y achievable at SLAC, it 1s possible to achieve much shorter wavelengths, 


including X rays. 


D. THE WAVE EQUATION 
Using the fields in 2-2 and assuming a slowly-varying amplitude and phase 
such that the terms with two derivatives may be ignored, we can write the optical 


wave equation as [3] 








i ai ar ~ 2 = + 1S (cos y, 0, 0) 
: c Ot 
2-15 


+ 2E B . 1 (-sin y, 0, 0) =-- "J, | 
c Ot " 


The total beam current Ji is the sum of all single-particle currents, 


J ,=-ecB 5°(x - 7), where the position of the ith electron is 7 and the transverse 


motion of the electron is determined 2-8. Separating the wave equation into its fast 


and slow varying factors and averaging over time gives 


A Pere n ) (e)—/ 
. i sega = 2V2REK (6) (8) Tessa, 2 
C 


Y 


The Bessel function factors describe the reduced coupling associated with fast 
periodic motion in each linearly-polarized undulator period. Following the field 
envelope traveling at speed c, the longitudinal coordinate becomes z = Z+ct, and 
simplifyies the differential operator from (0/dz + c' /6t) to c' G/Et. With respect 
to the light, the electrons then drift back with a relative speed c(1 - 8,). The 


summation in the wave equation can be replaced by an average < e > over 


sample electrons weighted by the electron particle density p. This leads to [3] 


Bp gn 2 22K UO AOE | oc re 


Ct y 


Finally, we introduce the dimensionless time t = ct/L, and the dimensionless 


current density j = 8N(enK(J,(&)-J,(4))L)’p/ymc, which simplifies the wave 


equation to 


da 


—= 7 (an me: 


Where the complex dimensionless field 1s a=4nNeK(J,(&) -J (§))LEe '*/2-y°me° and 


is consistent with the definition of the dimensionless field strength found below 
2-13. 

By assuming a plane wave, this derivation ignores the transverse evolution 
of the fields. When these effects are included, the three dimensional wave 
equation can be developed similar to the previous derivation and results in the 


parabolic wave equation [3] 


[-58 + S)aw RE v= - (yes) “ae 


x 


Here, V = & + oi and defines the transverse diffraction of the optical wave. The 


growth of the optical wave in the longitudinal direction 1s described by 2-19. 


E. ELECTRON BEAM QUALITY 
Of some concern to the FEL 1s the quality of the beam which feeds the 
undulator and determines its interaction with the light. One measure of the quality 
of the beam is the emittance. The emittance, €, is the average angular spread times 
the average beam radius, and the normalized emittance is €,= ey. The effect of 
electron beam emittance is to produce an effective velocity spread in the z 
direction which disturbs the FEL resonance condition [4]. As a minimum 


condition the transverse cross-section of the electron beam must be smaller than 
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that of the photon beam [5]. When we equate these conditions we arrive at 


ae 2-20 


where A 1s the shortest wavelength for which spatial coherence 1s desired [6]. 
Substitution of A= 4 nm shows a required normalized emittance less than €,=3.2 
mm-mrad. This is within the realm of current technology. If we project to the 
programs logical conclusion, achievement of hard X rays in the 0.15 nm regime, 


the required beam quality will be €, =| mm-mrad. 
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Il. SLAC X RAY FREE ELECTRON LASER 


A. SLAC PROPOSAL 

Live cell holography and high resolution imaging of DNA base pairs are 
two potential biological applications of an X ray laser [7]. The 2 to 4 nm 
wavelength arena corresponds to the "water window" where water becomes 
transparent to soft X rays. The duration of the pulse can be short enough that a 
biological sample will be unaffected by the radiation. 

The Free Electron Laser (FEL) appears to be the best possible source for an 
X ray laser. An advantage is that the FEL requires no light source and produces 
coherent radiation from noise in a process called Self Amplified Spontaneous 
Emission (SASE). Through this mechanism coherent X rays are created in a single 


pass which obviates the requirement for an X ray light source to feed the laser and 
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Figure 3-1. Schematic of the Slac linac driving a coherent X ray FEL. 
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the need for mirrors. Figure 3-1 is a schematic of the proposed system combining 
the SLAC linac with an undulator to produce coherent X rays. The system design 
would incorporate part of the SLAC linac to accelerate the electron beam to an 
energy of ymc’= 7 GeV with a peak current of I=2500 A and a bunch length of 
0.16 picoseconds. These pulses are injected into an undulator with a period of A, 
= 8 cm and a proposed length of 60 m, producing light of 4 nm wavelength and 
a peak power of 28 GW. Using three sections of the linear accelerator separated 
by two compression stages the electron bunch is accelerated from 10 MeV to the 


desired 7 GeV. Table 3-1 contains the design parameters for the X ray FEL 


<= 
= 
2 Se 
= aa 


| Beam Radius 0.0065 cm Peak Power 28 GW 










Table 3-1. Proposed parameters for an X ray FEL at SLAC 


We, 


planned at SLAC[1]. The peak current of 2500 A corresponds to a dimensionless 
current density / = 4000 and results in the FEL operating in the high gain regime. 
Most FEL's operate in the low gain regime where j~1. The evolution of the 
electrons through the undulator is governed by the pendulum equation 2-13, while 
the evolution of the optical amplitude and phase is determined by the optical wave 
equation 2-18 or 2-19. The dimensionless current determines the response of the 
optical wave to the electron bunching, coupling 2-13 and 2-18 or 2-19. Utilizing 
these equations, we initially examine the evolution of the electrons in phase-space 


(C, v). Figure 3-2 is a phase-space plot using the parameters of Table 3-1. At 


w** FEL Phase Space Evolution *** 
4324000 a,=0 Vom N=723 


Oe =3.6 


ln (1+P) 





Figure 3-2. Phase-space simulation where dimensionless current j=4000 
(I=2500 A). 
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startup (t=0), there is no optical field a = 0. The electrons are all started with 
random initial phases ¢ between -2/2 and 32/2 and periodic boundary conditions 
are assumed throughout the evolution.The initial electron phase velocities are 
spread randomly in a Gaussian distribution with standard deviation 6, = 3.6 and 
average v, = 0. Superimposed 1s a random exponential spread with characteristic 
width 6g = |. The final position of the electrons in the phase-space of Figure 3-2 
are bunched in the upper half of the separatrix at the end of the undulator (t=1). 
The electrons have overbunched through the undulator and are now taking energy back from 
the light. The change in phase velocity is related to the change in energy by Av 
= 4nNAy/y. Ideally, the electrons should be bunched in the lower half of the 
separatrix, indicating the electrons have given up their energy to the laser light. 
This is shown in the plot of the logarithmic power evolution from t=0 to | in the 
upper right hand corner of Figure 3-2. Again we see that the maximum power 
occurs at approximately t=*% and as the electrons begin to overbunch the power 
diminishes. As the electrons begin to bunch near € = 7/2, they drive the optical 
field phase 6, as shown by the graph in the lower right. This causes the separatrix 
to shift in phase space, so that the electron bunch is near the optimal relative phase 
C+o = x for gain. The simulations show that the laser saturates too early in the 
undulator and the optical power drops off before it can be utilized. In order to 
optimize the laser, we can decrease the length of the undulator. Figure 3-3 1s a 


phase space simulation with j=2000 which produces the same results as decreasing 


7 )| 


the undulator length by about 20%. Now we see that saturation occurs at the end 
of the undulator. The electrons are bunched in the lower half of the separatrix at 


the relative phase C+ =x where maximum gain occurs. 


*** FEL Phase Space Evolution *** 
4=2000 a=0 V=0 N=723 


G, = 3.6 G, = 10 





Figure 3-3. Phase-space simulation where the dimensionless current 
j=2000 ([=1250 A). 


Another solution to the early saturation and loss in power in the undulator 
is to taper the undulator by decreasing the undulator wavelength A, or decreasing 


the field strength B. To account for the force due to tapering, the pendulum 
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equation 2-13 takes the form [3] 


2 
ae la| cos(C +6) , oe 


4 
= 





dt 
where 6 is a constant torque given by 6 = -4nNK°AB/B(1+K") if we vary the 
undulator field strength and by 6 = -2xNAA,/A, if we vary the wavelength. Here, 


AB and Ad, are incremental changes in the undulator field strength and wavelength 


respectively. By setting €¢ = (ig k(t’) dt’ + kz - wt we account for the changing 
0) 


undulator wavenumber &,(t) = 27/A,(t). 

In Figure 3-4 the undulator simulated in Figure 3-2 with /=4000 is simulated 
using a taper of d6=SOnx. This corresponds to a AB/B = 3% change in the undulator 
field. Here the acceleration induced by tapering overcomes the natural deceleration 
of the untapered undulator which would extract optical energy accelerating the 
electrons. Comparing Figure 3-2 to Figure 3-4 we see that instead of reaching 
saturation three-fourths of the way through the undulator followed by a drop in 


power, now we continue to Increase power throughout the undulator. 
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Figure 3-4. Phase-space simulation with a tapered undulator. 


B. LONGITUDINAL FIELD DEVELOPMENT 

The optical field in the longitudinal direction can be generalized by 
following several sites along the complex wave envelope a(z). We scale our length 
in the z direction by the slippage distance NVA, the distance the electron pulse slips 
back relative to the optical pulse over one pass through the undulator. This 


corresponds to the distance over which the electrons and light exchange 
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information. The pendulum and wave equations 2-13 and 2-18 can be generalized 


to [3] 





d? 
- = |a.| cos (6... + 6) 
ate . 
3-2 
da | . 
=a Vane) 


As the electrons at each site pass through the undulator, they interact with a series 
of optical sites passing information from one site to the next. Figure 3-5 shows 
the effects of slippage and field growth. In the bottom left we see the exponential 
gain over the length of the undulator peaking at the end of the undulator. The 
bottom-middle illustrates the single-mode gain spectrum as a function of frequency 
around resonance. In the bottom-right window is the power development 
throughout the length of the undulator. Once light is produced, the power grows 
exponentially reaching its maximum at the end of the undulator. The middle-left 
depicts the growth in the amplitude of the light in a window that is two slippage 
distances in length -1 < z< 1. We see starting from noise, significant light 
develops near t = % through the undulator corresponding to the previously 
described power growth in the bottom-right. The middle and top-middle frames 
show the evolution of the power spectrum. Here we note the triangular tickmark 
indicates resonance v = 0. The rectangular tickmark indicates the wavelength of 


the final average power and shows development of power at wavelengths which 
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are higher than the fundamental. 


show the development of the electron phase velocity distribution. 
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Figure 3-5. Coherence evolution simulation showing the effects of slippage. 


decrease in the average electron phase velocity which corresponds to the 
development of high optical power. The final phase velocity spectrum (top-right) 
1s spread showing the acceleration of untrapped electrons to larger phase velocities 
and trapped electrons decreasing in energy. The pointed tickmark indicates the 


initial electron phase velocity at resonance while the rectangular tick mark indicates 
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Finally, the middle-right and top-right frames 


There is a 


the final average electron phase velocity showing a loss of energy. This shift in 
phase velocity accounts for the shift in the power spectrum showing optical power 


developing at longer wavelengths than the fundamental. 
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IV. THE FUTURE X RAY FEL AND ALTERNATIVES 


A. AFTER THE 4 nm FEL 

The 4 nm X ray laser as detailed in Chapter III may be feasible. We now 
look at the next step into the hard X ray regime. The follow-on experiment at 
SLAC has wavelength 0.15 nm and 1s detailed in Table 4-1 [8]. As discussed in 
Chapter II equation 2-20 the required normalized emittance (€,) to achieve full 
spatial coherence at A = 0.15 nm is 1 mm-mrad. Current technology can achieve 
3 mm-mrad at 5 GeV. The prognosis is good that RF guns with | mm-mrad at 25 


GeV will soon be available [8]. Assuming these parameters will be achievable 


ae i 
—— ae 
— a 


Table 4-1. Proposed parameters for 0.15 nm X ray FEL. 
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we now examine the viability of producing hard X rays in a FEL. Initial 
simulations show that saturation doesn't occur in the 70 m undulator outlined in 
Table 4-1. Instead, saturation occurs at 100 m as depicted in Figure 4-1 which is 
a phase-space simulation of a 150 m undulator. However, saturation yields peak 
power and for the purposes envisioned for the X ray FEL, power out 1s not a 


primary concern. So we need to examine whether coherent light 1s produced in 


*x*x*x FEL Phase Space Evolution *** 
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Figure 4-1. Phase space simulation for 150 m undulator utilizing the 
parameters outlined in Table 4-1. 
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a 70 m undulator. Figure 4-2 is a transverse multi-mode simulation examining the 
effects of optical guiding and mode distortion. Here, we see in the top-left the 
evolution of light along the length of the undulator and in the top-middle is a 
cross-section of the light at the end of the undulator. From this simulation we see 
that in an undulator of 70 m that while we do not achieve saturation we do 
produce light. The middle row of plots depict the position of the electrons in the 
beam. The bottom row depicts the average evolution in electron phase velocity 
which shows that at the end of the undulator the average phase velocity has shifted 
negatively indicating a loss of energy which has been given up to the production 
of light. But this simulation does not provide information as to whether the light 


is coherent along the longitudinal axis. 
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Figure 4-2. Transverse multi-mode simulation examining the effects of 
optical guiding and mode distortion. 
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In order to examine the coherence development in this design, again we 
look at a longitudinal multi-mode simulation in Figure 4-3. We find in the upper 
left column the development of a window of light through the undulator which ts 
representative of the development of light throughout the undulator. In the middle 
column we see that light 1s produced tn a narrow bandwidth indicating that the 
light 1s indeed coherent. The fractional spread 1n wavelength is AA/A = Av/2nN 
~ 9x10". In the bottom right corner is a plot of the dimensionless final power 


which equates to 170 MW of peak power. These results indicate that upon 
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Figure 4-3. Longitudinal multi-mode simulation depicting coherence 
development. 
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production of an RF gun which meets the requirements of Table 4-1, the move to 


hard X rays can be made. 


B. AN ALTERNATIVE APPROACH. 

The design for a 4 nm X ray FEL outlined in Chapter III appears to be 
viable. However, undulator technology is expensive and the requirement for a 60 
m undulator could be quit costly. If the undulator could be shortened and still 
achieve the same results, that would be a prudent move. An option for boosting 
the gain in an undulator where we are working with weak optical fields 1s to add 
a dispersive section producing a klystron. Here we explore the production of soft 


X rays with A = 4 nm in a klystron. 


——— > a 


a 
a ee 
0.0065 cm Wavelength A 0.15 nm | 


Table 4-2. Proposed parameters for a 4 nm klystron FEL. 
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In the first undulator section of the klystron, called the "modulator", the 
electrons are prepared for bunching. While they traverse the dispersive section 
without interaction with the optical fields, the electrons become bunched. As the 
enter the second undulator, called the "radiator", the bunched electrons begin to 
radiate producing coherent light. Table 4-2 lists the parameters for this proposed 
design. The strength of the dispersive section 1s described by the dimensionless 
value D [9]. The effect of the dispersive section is described by AC = vD. The 
dispersive section strength D was determined by a series of phase space 


simulations which revealed the most gain was produced with D=0.5. Again to 


wkekke PEL Wrap Evolution **t#*# 
j=932 V5=0 a=0 
D#=0.5 N=413 Og=1 .04 Og=0. 26 


9MMMM71. Pv, 2) 


fa(z,T) | 0 


cy eel ingisey TBS 





Figure 4-4. Simulation depicting coherence development for klystron X 
ray laser. 
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determine the coherence of the light produced we inspect the results of the 
longitudinal multi-mode simulation in Figures 4-4 and 4-5. Of primary concern 
is the middle window in both figures which shows a spike of very narrow 
bandwidth indicating good coherence development. Both figures are displayed in 
order to demonstrate that by changing the random number generation in the 
simulation we randomize the position of the electrons at startup. While the final 
shape 1s somewhat altered the development of coherence is still observable. These 
same variations have been tested with the proposed SLAC X ray laser producing 


the same coherence results. The right column displays the average electron 
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Figure 4-5. Simulation of klystron X ray laser varying electron starting 
position. 
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evolution showing maximum loss of energy of the electrons at the end of the 
radiator. Finally, the power presented in the bottom right equates to 2 GW of peak 
power. By comparison to the proposed design in Chapter III, we can produce 
coherent X rays with a 17 m undulator instead of a 60 m undulator which would 


seem to be advantage in reducing expenses. 
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V. CONCLUSIONS 


The theory for the FEL which has been presented applies not just to the X 
ray FEL specifically but to any FEL in general. The tunability which has been 
shown makes the FEL suitable for many military and industrial applications. The 
FEL as a shipboard, land based or airborne weapon can respond to a threat at the 
speed of light and does not require replenishment of a finite magazine. The FEL 
will have a significant impact upon the military of the future and the defense of the 
nation. 

The proposed X ray FEL 1s the first viable entry into the regime of coherent 
X ray light sources. Using the proposed SLAC FEL parameters in Table 3-1, 
simulation results indicate that X rays in the 4 nm regime are achievable now. The 
future regime of 4 = 0.15 nm requires an electron beam energy of 25 GeV, a peak 
current of 5000 A, and A, = 4 cm which are achievable. However, the required 
beam quality of €,= 1 mm-mrad 1s not available, but is expected in the near future. 
This opens the door into a new exciting biological arena for the Free Electron 
Laser. In Table 5-1, we compare our results with data generated in a coherence and 
linewidth study by Lawrence Berkeley Laboratory and Lawrence Livermore 
National Laboratories [10]. The number of undulator periods necessary to reach 
saturation, N,,,, 1s plotted as a function of the dimensionless current density j. In 
this study, LBL and LLNL used the GINGER particle simulation code to evaluate 
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the effect of current on the FEL. We see that the simulations we report are 


consistent with their findings. 


XLBL/LLNL |! NPS 


400 


200 
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Table 5-1. Data comparison of LBL/LLNL and NPS simulations. 


Finally, as an alternative to the currently proposed SLAC X ray laser 
utilizing a 60 m undulator, we have presented a variation utilizing a klystron. By 
doing so, the required undulator length is reduced to only 17 m. We have 
presented simulation results which makes this a feasible and potential cost savings 


option which requires further investigation. 
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